Motion of a heavy impurity through a Bose-Einstein condensate. 
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Abstract 

We study motion of a point-like impurity in a Bose-Einstein condensate at 

T = 0. By solving the Gross-Pitaevskii equation in a perturbative manner 

we calculate the induced mass of the impurity and the drag force on the 

impurity in 3D, 2D and ID cases. The relationship between the induced 

mass and the normal mass of fluid is found and coincides with the result of 

the Bogoliubov theory. The drag force appears for supersonic motion of the 

impurity. In ID the drag force is investigated also on the basis of the exact 

Lieb-Liniger theory, using the dynamic form factor, which has been evaluated 

by the Haldane method of the calculation of correlation functions. In this 

theory the force appears for an arbitrarily small velocity of the impurity. The 

possibility of measuring the form factor in existing experiments is noted. 
PACS numbers: 03.75.Kk, 03.75.Lm, 67.40.Yv 
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I. INTRODUCTION 



One of the most important peculiarities of Landau theory superfluidity is the existence 
of a finite critical velocity. If a body moves in a superfluid at T = with velocity V 
less then Vc, the motion is dissipationless. At V > Vc a drag force arises because of the 
possibility of emission of elementary excitations. However, both theoretical and experimental 
investigation in superfluid ^ife are difficult. The critical velocity in ^ife is related to creation 
of rotons, for which one has no simple theoretical description. Further, an important role is 
played by complicated processes involving vortex rings production. 

The situation in low-density weakly-interacting Bose-Einstein condensed (BEC) gases is 
simpler. The Landau critical velocity in this case is due to Cherenkov emission of phonons 
which can be described by mean-field theory. Due to the presence in the theory of an 
intrinsic length parameter - the correlation length ^- the friction force for a small body does 
not depend on its structure. Vortex rings in the BEC cannot have radius less then ^ (see 
[1]) and it is reasonable to believe that probability of their creation by a small body is small. 
Thus quantitative investigation of critical velocities in BEC are very interesting and can be 
used to probe the superfluidity of a quantum gas. 

Recently existence of the critical velocity in a Bose-Einstein Condensed gas was confirmed 
in a few experiments. At MIT a trapped condensate was stirred by a blue detuned laser 
beam [2] and the energy of dissipation was measured. The critical velocity was found to be 
smaller than the speed of sound due to emission of vortices. The diameter of the laser spot 
in this experiment was of a macroscopic size and was large compared to the healing length. 
An improved technique allowed measurement of the drag force acting on the condensate in 
a subsequent experiment [3]. 

The analytical study of flow of the condensate over an impurity is highly nontrivial due 
to the intrinsic nonlinearity of the problem arising from the interaction of the particles in the 
condensate. In one dimension the dissipation could occur at velocities smaller than predicted 
by Landau's approach due to emission of solitons [4] . The dependence of the critical velocity 
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on the type of the potential was studied both by using a perturbative approach and numerical 
integration in [5,6]. The effective two dimensional problem was considered in [7]. In this 
work generation of excitations in the oscillating condensate in a time dependent parabolic 
trap in the presence of a static impurity was studied analytically. A three-dimensional 
flow of a condensate around an obstacle was calculated numerically by integration of the 
Gross-Pitaevskii (GP) equation and emission of vortices was observed [8,9]. 

In this letter we study a 5-function perturbation moving at a constant velocity in a 
condensate. We find analytically the depletion of the superfiuid fraction and the drag force. 

II. THREE-DIMENSIONAL SYSTEM 

Let us consider an impurity moving through a three-dimensional condensate at T = 0. 
One of the possible reahzations of this model could be scattering of heavy neutral molecules 
by the condensate. 

A. Effective mass and normal fraction 

We start from the three-dimensional energy functional of a homogeneous weakly- 
interacting Bose gas in the presence of a 5-function perturbation (an impurity) moving 
with a constant velocity V 

where -0 is the condensate wavefunction, fj, is the chemical potential, m mass of a particle 
in the condensate, g — Anh'^a/m and Qi — 27r^^6/m are particle-particle and particle- 
impurity coupling constants, with a and b being the respective scattering lengths [10]. We 
will assume that the interaction with impurity is small and we will use perturbation theory. 
By splitting the wave function into a sum of the unperturbed solution and a small correction 
V^(r, t) = 00 + ^'0(r, and linearizing the time-dependent GP equation with respect to 6ip, 
we obtain an equation describing the time evolution of Sip 
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+gi5{r-\t)(Po. (2) 

In a homogeneous system 0o is a constant fixed by the particle density 4>o — ^/n and 
II — gn — mc^. 

The perturbation follows the moving impurity, i.e. Sif) is a function of (r — Vt), so the 
coordinate derivative is related to the time derivative ddip{r — Yt)/dt — —YVSiplr — Yt). 
We shall work in the frame moving with the impurity r' = r — and the subscript over r 
will be dropped. 

Eq. (2) for a perturbation in a homogeneous system can be conveniently solved in 
momentum space. In order to do this we introduce the Fourier transform of the wave 
function Sijj]^ = J e~'^^'''5ijj{r) d^x. Eq. (2) becomes 

(hk-\-^-mc^^ SiP^-mc''{5iP_y,y-gi<l>o^O. (3) 

Substitution of A; ^ —k and complex conjugation of (3) give the second equation. The 
obtained system of linear equations can be easily solved 



Let us calculate the energy of the perturbation. By neglecting terms of the order of 

gilSi/jl"^ the energy functional (1) becomes 

E^Eo+g^4+ lhk.\\S^lj^\'-^^ + ^{5ij+5r)r=o. (5) 

Here Eq — Ngn/2 is the energy of the system in absence of the perturbation. We 
expand the 3^'^ and 4*^ terms in powers of V. To avoid the largc-fc divergency, one must 
also introduce a renormalization of the scattering amplitude. It is sufficient to express the 
couphng constant in the 2"^^ term of eq. (5) in terms of the scattering amplitude b using the 
second order Born approximation: 



^^ = ^^ + ^7 [-2^) ■ 



Finally, by carrying out the integration over momentum space and considering A^j^p 
impurities with a concentration given by % = N^^p/N we obtain the energy per particle 
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If we set = we recover Bogoliubov's corrections to the energy in the presence of 
quenched impurities [11,12]. Note that even if the "mean-field" energy obtained from the 
GP equation in the absence of impurities (x = 0) leaves out terms of the order of (na^)^/^, 
the equations we obtain in the presence of impurities in a perturbative manner still correctly 
describe the effect of the disorder up to the terms of the order of (na^)^/^. 

If y 7^ a quadratic term in the impurity contribution to the energy is present. It can 
be denoted as xm*V^/2 with 
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m* = h^(na'y/^-] m (8) 



being the induced mass, i. e. the mass of particles dragged by an impurity [13]. Applicability 
of the perturbation theory demands m* to be small compared to m. This gives the condition 
{na^y^^ <^ 1. At zero temperature the interaction between particles does not lead 
to depletion of the superfiuid density and the suppression of the superfluidity comes only 
from the interaction of particles with impurities. Thus (8) defines the normal density p„ = 
Xpm*/m. 

This result is in agreement with the one obtained by the means of Bogoliubov transfor- 
mation starting from the Hamiltonian written in the second-quantized form in the presence 
of disorder [11,12]. 

The normal density of a superfluid is an observable quantity. It was evaluated in liquid 
"^He by measuring of the moment of inertia of a rotating liquid or by measuring of the second 
sound velocity. Both methods can be, in principle, developed for BEG gases. 
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B. Drag force and energy dissipation 



The force with which the impurity acts on the system is 

F = -/ |V'(r)rVM(r))d^x = ^,(V|V'(r)r)r=o- (9) 

Expanding the wave function into the sum of 0o ^-nd Sip and neglecting terms of order 
Sip"^ we obtain 

ik [SilJi, + (dip-k)*] 

2{gicj>ofik(h^e/2m) 



~ J {hk-V + iO)2 - ^ + 2mc2) (27r)' ' 

where we added an infinitesimal positive imaginary part +iO to the frequency k • V according 
to the usual Landau causality rule. The drag force is obviously directed along to the velocity 
V. 

We can carry out the integration with respect to cosi?, where i? is the angle between the 
momentum k and velocity V, using the formula -^^^ inS^x). Due to the integration 

between symmetric limits, only the imaginary part contributes to the integral with respect to 
cos '&. The poles in the integration over cos i? appear if the square root in the denominator is 
smaller than one, which leads to the restriction on the values of momentum which contribute 

< k^^ = 2m{V^ - c^f'^/n . (11) 

Thus the energy dissipation takes place only if the impurity moves with a speed larger than 
the speed of sound. Integration with respect to A;, taking into account restriction (11), finally 
gives 

Fv = 47rn6'my2(l - /V^f . (12) 

The energy dissipation, E — —FyV, can be evaluated by measuring the heating of the gas. 

For large V the force is proportional to V^. The energy dissipation per unit time can 
then be presented &s E — —jE with the damping rate 7 ~ nb^V. 
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Note in conclusion that our perturbative calculations can not describe processes involving 
dissipation of energy due to creation of quantized vortex rings. Such a creation is possible 
at V < c but has a small probability for low velocity and for a weak point-like impurity. 

III. LOW DIMENSIONAL SYSTEMS 

In this type of experiment the role of the impurity can also be played by a laser beam 
with small enough size and intensity. The Fourier components of the perturbed wavefunction 
^■^k are given by the formula (4), which is derived in an arbitrary number of dimensions. 
The only difference is in the substitution of d^k/{27r)^ with d^/{27r)^ in the integrals. 

A. Two-dimensional system 

There are different possible geometries of the experiment. One can create a two- 
dimensional perturbation in the three-dimensional condensate. Such a two-dimensional im- 
purity can be created, analogously to the MIT experiment [2,3], by a thin laser beam. Such 
a beam creates a cylindrical hole in the condensate, which is stirred by moving the position 
of the laser beam. Another possibility is to fix the position of the laser beam along the long 
axis of an elongated condensate, so that the dissipation can be studied by shaking the trap 
and exciting the breathing modes. The problem is to create a beam with a diameter which 
is small with compared to the correlation length. The theory can be easily generahzed for 
beams of finite diameter. The intensity of the beam can be tuned to satisfy the condition 
of a weak perturbation. 

The more interesting possibility is the investigation of true two-dimensional condensates, 
which can be created in plane optical traps, produced by a standing fight wave. If the light 
intensity is large enough, tunneling between planes is small and the condensates behave as 
independent two dimensional systems. The impurity can again be created by a laser beam 
perpendicular to the condensate plane. Another possibility is to use impurity atoms, which 
can be drive by a laser beam, with a frequency close to the atomic resonance of the impurity. 
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In two dimensions one has (Pk — kdkd'd. The drag force F is different from zero only if 
the denominator has poles, which means that the velocity V must be larger than the speed 
of sound c. Only momenta smaller than kmax (see eq. (11)) contribute to the integral. 

Integration for the velocities V > c gives us the drag force in the two-dimensional case: 

Fr^g^n,m\V'-^)/{n'V), (13) 

where n2 is two-dimensional density. 

In a quasi two-dimensional system, i.e. when the gas is confined in the z-direction by the 
harmonic potential (1/2) mw^r^, the two-dimensional coupling constant equals gi = \/27r^^, 
where az — ^Jh/mujz is the oscillator length and h is the three-dimensional scattering 
length. (We consider here only the mean- field 2D situation. See [14], § 17 for a more detail 
discussion.) 

The calculation of the effective mass gives 

m* = g1n2m/{A-K'^h'^c^) . (14) 

Notice again that our calculations do not take into account creation of vortex pairs which 
is possible at y < c. 



B. One-dimensional system. Mean-field theory 

In one dimension the integration is straightforward. The integration (10) over k gives 
2m if y > c and zero otherwise. So, the force is F^P — 2gfnim/h'^ , where ni is the hnear 
density. In a quasi one dimensional system (i.e. a very elongated trap or a waveguide) there 
are no excitations in the radial harmonic confinement and the coupling constant is given by 
gi — — ^^^^ with hiD — —a]_/b, where a± — \Jh/mcu± and one has F^^ — 2nih^ /mbfjj. 

An interesting peculiarity is that the result does not depend on the velocity V (where, 
of course, the velocity must be larger than the speed of sound). The same result for the 
5-potential was found in [6]. Calculation of the effective mass gives m* — g^ni/2hc^. 
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In a ID system energy dissipation is possible ai V < c due to creation of the "gray 
solitons" first considered in [15]. Non-Hnear calculations [4] show that the critical velocity 
for this process decreases with increasing coupling constant gi. 

This theory can be checked in an experiment in a three-dimensional condensate. The 
two-dimensional impurity can be presented by a moving light sheet. 

C. One-dimensional system. Bethe-ansatz theory 

We saw in the previous subsection that for a weakly interacting impurity the drag force 
appears only when the impurity velocity V is larger that the Landau critical velocity, which 
is equal to the velocity of sound c. The situation is, however, different in the Bethe-ansatz 
Lieb-Liniger theory of a ID Bose gas [16]. According to this theory excitations in the system 
actually have a fermionic nature. Even a low frequency perturbation can create a particle- 
hole pair with a total momentum near 2pp = 2fikF = h27Tni. To calculate the drag force for 
this case we will use the dynamic form factor of the system a {u!,k) (we follow notation of 
[17], § 87). The dissipated energy at T = can be calculated as 

-oo 

where U{u!, k) — 2'Kgi5{u! — kV) is the Fourier transform of the impurity potential U{t, z) — 
giS{z — Vt). One has \U (a;, k)f = 27:gftS{u! — kV), where t is "time of observation". Thus 
the energy dissipation per unit of time is 

E = -FyV = J^ka {kV, k) , (16) 



where Fy is the drag force. We will try to estimate the velocity dependence of Fy- 

For low frequency dissipation the important values of k are near 2kF- According to [18] 

a(c;,2A;F) ~c^(''-2),cu ^ 0, (17) 

where ri = = ^^^"^ > 2 is the characteristic parameter of a ID Bose gas. In the mean- 
field limit when ni — > oo the parameter 77 ^ oo. In the opposite case of a small density 
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bosons behave as impenetrable particles (Girardeau limit [19]) and the dynamic form- factor 
coincides with the one of an ideal Fermi gas. In this limit rj — 2. 

In the general case one can calculate a {u!,k) at small cu and k ~ 2A;i7 generalizing the 
method of Haldane [20] for the case of time-dependent correlation functions. Calculations 
give 



nic f hu \ ^ J. ^cA/c' 



^i^,k) = -^[^^^j f[-—] ,u;>0,k>0, (18) 

where k — 2k f + Ak and the function f{x) is 

f{x)^A{n){l-x'y^"~" (19) 

in the interval \x\ < 1 and is equal to zero at \x\ > 1 (see also [21]). The constant A{ri) 
can be calculated in two limiting cases: A{ri = 2) = ix/A (see [14] § 17.3) and A{ri) ^ 
SttV [(8C)T2 (i)] , where C = 1.78... is the Euler 's constant, for 77 > 1 (details of the 
calculation will be published elsewhere). 

Substituting (18) into (16) we finally find velocity dependence of the drag force: 



20Fr(2±i) hv V'c 

Equation (20) is valid for the condition V ^ c. 

Thus in the Girardeau strong-interaction limit Fy ^ V and Bose gas behaves, from the 
point of view of friction, as a normal system, where the drag force is proportional to the 
velocity. On the contrary, in the mean-field limit the force is very small and the behavior of 
the system is analogous to a 3D superfluid. However, even in this limit the presence of the 
small force makes a great difference. Let us imagine that our system is twisted into a ring, 
and that the impurity rotates around the ring with a small angular velocity. If the system is 
superfluid in the usual sense of the word, the superfluid part must stay at rest. Presence of 
the drag force means that equilibrium will be reached only when the gas as a whole rotates 
with the same angular velocity. From this point of view the superfluid part of the ID Bose 
gas is equal to zero even at T = 0. Notice that in an earlier paper [22] the author concluded 
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that ps = p at T = for arbitrary 77. We believe that this difference results from different 
definitions of Ps and reflects the non-standard nature of the system. 

Equation (20) is equivalent to a result which was obtained by a different method in [23], 
with a model consisting of an impurity considered as a Josephson junction. Notice that the 
process of dissipation, which in the language of fermionic excitations can be described as 
creation of a particle-hole pair, corresponds in the mean-field limit to creation of a phonon 
and a small-energy soliton. It seems that such a process cannot be described in the mean- 
field approach in the linear approximation. 

Experimental confirmation of these quite non-trivial predictions demands a true one- 
dimensional condensate, where non mean-field effects can be sufficiently large. Such con- 
densates have been investigated for the first time in experiments [24,25]. In experiments 
[26,27] condensates have been created in the form of elongated independent "needles" in 
optical traps, consisting of two perpendicular standing laser waves. The role of an impurity 
in this case must be played by a light sheet, perpendicular to the axis of condensates and 
moving along them. 

Notice also, that application of the additional light waves in this experiments of this type 
allows one to create a harmonic perturbation of the form 

U{t, z) = Uo cos {ujt -kz),k^ 2kF + M (21) 

with small cu and Ak. Such potential with was used in [26,27] for experiments with ID 

condensate in a periodic lattice. However, for a small amplitude Uq, measurement of the 
dissipation energy Q gives, according to (15) , the dynamic form-factor S {cu, k) directly. 

IV. CONCLUSIONS 

We have studied motion of an impurity through the condensate at zero temperature by 
considering the perturbation of a stationary solution of the GP equation. We calculated the 
induced mass which contributes to the mass of normal component. We find that the motion 
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at small velocities is dissipationless in one-, two-, and three- dimensional systems, although 
movement with velocities larger than the speed of sound leads to a non-zero drag force due to 
Cherenkov radiation of phonons. The expressions for the drag force are calculated. We used 
results for the dynamic form factor of exact Lieb-Liniger theory to investigate the velocity 
dependence of the drag force in a ID system. The form factor was calculated with the 
help of the Haldane method of calculations of correlation functions. The drag force exists 
at an arbitrarily small velocity of motion, but is very small in the mean-field limit. The 
dynamic form-factor can be also directly measured by applying a harmonic time-dependent 
perturbation on one- dimensional condensates [26,27]. 

Acknowledgments: We thank T. Esslinger, F. D. M. Haldane, B. Jackson, Yu. Kagan, 
V. Korepin, E. Sonin and S. Stringari for useful discussions. 
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